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Abstract: It is theoretically and experimentally shown that photons 

emitted by statistically independent incoherent classical light sources and 
measured in the far field in spatially separated modes may display spatial 
correlations akin to path-entanglement of photons produced by quantum 
sources. By measuring higher order photon-correlations at different loca¬ 
tions, i.e., m photons in one mode and one photon in another mode, we 
experimentally demonstrate for m> 6 a violation of Bell-type inequalities 
for spatial degrees of freedom. The spatial correlations among the photons 
can be understood from state projection where the detection of the first m 
photons projects the sources onto a state which emits the subsequent photon 
in a strongly correlated manner. From this perspective the entanglement and 
violation of Bell’s inequalities appears as a consequence of nonvanishing 
cross correlations between noncommuting quadrature phase components of 
the two spatially separated fields after m photons have been recorded. In this 
way we show that classical systems may produce spatial field correlations 
violating local realistic theories. 
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1. Introduction 

Since Bohr formulated his correspondence principle [1] it is arguably agreed that classical 
physics can be used to explain quantum phenomena in the classical limit, i.e., in the case of 
high quantum numbers. For example, electron wave functions reduce to Bohr orbits in the limit 
of high principal and orbital quantum numbers, and laser amplitudes display vanishing fluctua¬ 
tions in the limit of an infinite mean photon flux. By contrast, it is less anticipated that quantum 
mechanical predictions, i.e., implications considered to be purely quantum mechanical, can be 
fruitfully applied for a better understanding of classical systems. In this context it was recently 
proposed that a violation of Bell’s inequalities, so far considered to be a typical quantum me¬ 
chanical phenomenon displayed only by quantum systems, may be used to quantify correlations 
between coupled degrees of freedom of classical systems [2, 3, 4, 5, 6, 7, 8, 9, 10]. Such strong 
correlations of classical systems, called classical entanglement or nonquantum entanglement 
in the literature [5, 8], have been used to describe the specific coherence properties of classical 
light beams, either for discrete or continuous variables [1 1, 5, 7, 8, 9]. However, classical entan¬ 
glement has been demonstrated so far only for coupled degrees of freedom of the same particle, 
e.g., amplitude and polarization of a photon in a classical light beam, and was expected to reach 
its limits if true nonlocal multiparticle entanglement is considered. It has thus been conjectured 
that nonlocal multiparticle entanglement is of exclusive quantum nature [2]. In this letter we 
demonstrate that Bell’s inequalities can be equally violated with spatially separated particles 
produced by classical systems, i.e., photons emitted by classical light sources recorded in spa¬ 
tially separated modes. We thus show that nonlocal multiparticle entanglement is not restricted 
to the quantum world but may be observed also with classical systems. 

The paper is organized as follows: In Sect. 2 we present the system under investigation and 
introduce the normalized first and second order spatial intensity correlation functions employed 
for Bell’s inequalities involving spatial degrees of freedom. In Sect. 3 we expand the analysis 
and introduce intensity correlation functions of higher order which we use, with a new inter¬ 
pretation, to violate the Bell inequalities. In Sect. 4 we present our experimental results and in 
Sect. 5 we finally conclude. 

2. System under investigation 

Recently, we demonstrated that photons emitted by statistically independent quantum sources, 
e.g., single photon emitters (SPE), may be path-entangled if the photons are recorded in the 
far field such that the individual photon source remains unknown [12]. The entanglement was 
proven by a violation of Bell’s inequalities [13, 14] formulated for spatial degrees of freedom. 
The question remains whether such nonlocal correlations among photons recorded in spatially 
separated modes can be produced with classical light, i.e., photons emitted by classical sources. 
To investigate this question we generalize in what follows the system of two SPE discussed in 
[12] to classical light sources. 

We thus consider a system of two statistically independent identical classical light sources, 
e.g., thermal light sourcess (TLS), located at R/, I = 1,2, and two detectors placed in the far- 
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Fig. 1 . Scheme of the investigated setup. Two statistically independent Thermal Light 
Sources (TLS) at Ri and R2, separated by a distance d, emit photons which are registered 
by two detectors Di and D2 in the far-field at ri and r2, respectively. Assuming that each 
detector registers coincidentally exactly one photon the second order correlation function 
’ *”2) measured. 


field of the sources at r^, j = 1,2 (see Fig. 1). Due to the far field condition and the indistin- 
guishability of the emitted photons the electric field operator (r^) is written in the form 

£W(ry) = ^0 (fli +e'^id2) = (E^^Hrj)y , (1) 

where (So is the amplitude of the total electric field at and 5j = 5j{rj) = = kd sin(0j) 

describes the relative optical phase accumulated by a photon emitted from the source at R 2 
with respect to a photon emitted from the source at Ri, both reaching the y'-th detector at (cf. 
Fig. 1 ). Here, A: = ^ is the wave number of the photons, the operator d; denotes the annihilation 
of a photon from source I and the electric field is described by a scalar since we only consider 
a single polarization. 

Starting from the Bell inequalities for continuous variables [14] 

— XY<xy — xy'+x'y+x'y'—x'Y—yX<0, (2) 

which hold for all x, x' <X and y, y' < Y, we can, by setting X = Y = 1, interpret .x,x',y,y 
as probabilities and use Glauber’s first and second order intensity correlation functions [15] as 
input to Eq. (2). Hitherto the spatial intensity correlation functions of hrst and second order at 
equal times, defined as [15] 

G(')(ri)=(£-(ri)£+(ri)) 

G(^^(ri,r2) =(£^(ri)£'^(r2)£+(r2)£+(ri)) , 

have to be cast into a form as to be interpreted as detection probabilities of single photon detec¬ 
tion events 0 = G^'^ (ri) and joint two-photon detection events (ri, r2)ocG(^)(ri,r2), 

respectively [12]. 

If both light sources are incoherent and statistically independent the emitted intensity is dis¬ 
tributed isotropically in space. The first order correlation function G^'^rj) is thus a constant. 






independent of ri. By contrast, the second order correlation function for two statistically inde¬ 
pendent SPE reads [16] 

=2 ^0 (1-fcos(5i-52)) , (4) 

whereas for two statistically independent TLS with identical mean photon numbers = («) 
it is given by [16] 

^2^7’Ls(*’i>*‘2)=6^0^(n)^^l + ^cos(5i-52)^ . (5) 

In order to work with a probabilistic theory that remains valid (i.e., normalized) for arbitrary 
light sources, we introduce in what follows two more detectors in the setup where the first pair 
of detectors is placed at arbitrary positions 5i and 52 and the second pair at positions ni=5i + n 
and 7^2 = 52 -f With the four detectors at 5i, 52,7ti and 7^2, it is possible to measure altogether 
six second order correlation functions. 

(21 

In case of two SPE the six second order correlation functions read (cf. Eq. (4)) 

= 2 (1 -f 5 />£COs(5i — 82 )) 

^2^5/>£(^1’^2) = 2 (1 -f SPE COs{5l — 82 )) 

^ 2 ^ 5 />£(^i’ ^ 2 ) = 2 S’Q{l — 'f '2 specos{ 8 i — 82 )) (6) 

= 2 S’Q{l — y2 SPECOs{8i — 82)) 

G2^5/>£(5i, TTi) = G^'Ipp { 82 , 712 ) =2S'Q{l — 'f2 SPe) , 
where the visibility ^ spe has been added artihcially to incorporate possible experimental im¬ 
perfections [12]. By summing over all six correlation functions we obtain the following nor¬ 
malization factor 

L G^ilpEi^i^ ^j) + G^?spe { 81 , 82 ) + G^^Ipe (^ 1 , ^ 2 ) = 2 £’^{6 - 2^2 spe) , (7) 

U=1.2 
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^2 SPE (81,7^2) 

^2 

^2 SPE {81,7^1 ) —£’2SPe{^^^ 

whereas the probabilities for single photon detect 
rived from the two-photon detection probabilities 

p(l) (-x. W pP) 


-[i+r2S/.£COS[di-d2)) 

O — Z ^2 SPE 

7- 7^ -(1 -'^2S/>£COs(5i -52)) 

0 — 172 SPE 

■(1 — >2 speoos{8\ — 82)) 

1 — SPE 


P2 SPE {81) — 
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P2 SPEi8l,(7(') 

a= 

= 7rj ,712,52 


Pi'Li82) = 

E 

P2 SPEi^’82) 


6 — 2'f2 SPE 

p(2) ^ s „ ^ ^-'^2 SPE 

P2SPe{82,7Z2)- , 

rvtrvr. detectioo events P^‘jp£(5i) andP2^‘^^£(52) are de- 
lilities via 

3-^2 SPE _ 1_ 

6 — 2% SPE 2 


a=K\ 


3-^2 SPE _ 1_ 
6 — 2^2 SPE 2 


(9) 
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Identifying in Eq. (2) x = x' = ^2 y = / = PiIpe^^i) 

the probability to measure a single photon at position 5i, 5[, &2 and 5!^, respectively, and the 

CT) 

products xy,... as the joint two-photon detection probabilities we arrive at 

the following set of Bell inequalities 


-1< . (cOs(gi-g2)-COs(gi-g^) 

0 — 172 SPE ^ 

-l-cos(5i' — ^ 2 ) -l-cos(5i' — 52)1 -I- - —^ 1 < 0 . 

/ 0 — 172 SPE 


( 10 ) 


Using for the arguments of the cosine-functions in Eq. (10) the upper bound Bell angles 7r/4, 
37r/4,7r/4,7r/4 [17] for testing the upper bound, and the lower bound angles 37r/4,7r/4, 37r/4, 
37r/4 for testing the lower bound, one obtains a violation of Eq. (10) for different visibilities 
Pi SPE- Since the term in brackets obtains maximally (minimally) the value 2v/2 (—2^/2) for 
the upper bound angles (lower bound angles) the upper bound can be violated for spe > 
2/(1 + v/2) ~ 0.83 and the lower bound for ^2 spe > 1 /v^ ~ 0.71. Note that for ^2 spe = 1 

( 2 ) 

we obtain from Eqs. (8) and (9) for the two-photon detection probabilities Pj spe 


PPsPEiSu^) 

= ^(1+cos(5i-52)) 

Pi 

= ^(1+cos(5i-52)) 

^2 ^2) 

= ^(1-cos(5i-52)) 

Pi SPE^^^i^) 

= ^(1-cos(5i-52)) 

Pi 5/>£(^li ^1) 

“ Pi SPE^^’^^) = 0 j 


and for the single photon detection probabilities 



L 

Pi SPEi^i, Ct) — 

a= 

= 7Ci ,7r2,52 


P^%EiSl) = 

E 

Pi spe(^> = 


a=7i:i ,7r2,5i 

so that we are left with the more familiar Bell inequality [14] 


( 11 ) 


( 12 ) 


— 1 < ^ ^cos(5i — 52) — cos(5i — 52) -l-cos(5( — 52) -l-cos(5{ — 5^)^ — ^ < 0 . 


(13) 


In what follows we employ the same approach as above to study whether a violation of the 
Bell inequalities Eq. (2) can be obtained also for photons emitted by two statistically indepen¬ 
dent TLS. To this goal we derive by use of Eq. (5) the following set of second order correlation 
functions 


r(2) 

'-'1 TLS 

(5i,52) 

'-0 

II 

(«; 

>^{1+^1 

TLS cos (5i — 

52)) 

^(2) 

'^1 TLS 
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52)) 

^(2) 
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II 
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i^(l-r2 

T LS COS {5i — 

52)) 

r(2) 

'-'1 TLS 


“ P'1 rz,s(^> ^ 2 ) — 

6 5]/(«)"(! 

-Pi 


(14) 





where, in analogy with Eq. (6), we introduced again artificially the visibility tls to incorpo¬ 
rate experimental imperfections and also to infer the least visibility required to violate the Bell 
inequalities. Again, the sum over all second order correlation functions yields the normalization 
factor transforming the correlation functions into probabilities, which this time calculates to 

E Gf\^s(5i,nj) + Gfj^^si^u52) + Gfj^{ni,n2) = 6 4{nf{6-2y2TLs)- (15) 


and 


f2^rz.s(^i>^2) (l-l-f2rL.soos(^i-^)) 

O — 2 >2 TLS 

(1 - TLS cos (5l - 52 )) 

O — 2 >2 TLS 

-(1 — 7-L5COs(5i — 52 )) 


=2 —-(1 “ rL5COs(5l - 

b — 172 TLS 

(50 = (52) = = - . 

2TLS\ w ITLSy 2) 2 

from Eq. (16) and P2'^^^(50, P^tls^^) from Eq. (17) 
^cos(5i — 52 ) — cos(5i — 50 


Plugging pfL(5i,52) 
obtain 

_Y ^ TLS 

“ 6 — 2^2 TLS 

+ cos(5[ — 52 ) +cos( 


into 


(17) 
Eq. (2) we 

(18) 


Since Eq. (18) is structurally identical to Eq. (10), the upper and lower bound can be be vio¬ 
lated for identical visibilities as in Eq. (10), i.e., for tls > 2/(1 -f s/l) « 0.83 for the upper 
bound and lO tls > 1 /v/2 ~ 0.71 for the lower bound. However, the visibility of the second 
order correlation function for two photons emitted from two statistically independent classical 
sources acquires at most a value 'P =\/ 2 , which is the case for two statistically independent 
coherent sources =1/3 for two statistically independent TLS), even in the ideal case of 
negligible experimental insufficiencies [18]. It is obvious that with this value for the visibilities 
the Bell inequalities of Eq. (18) can not be violated. 


3. Violation of Bell’s inequalities using higher order intensity correlations 

In this section we investigate whether by recording m -f 1 photons (m > 2) emitted from two 
statistically independent TLS we can produce a visibility 'P^'tls'^ enough to violate Eq. (2). 
Hereby, we assign via post-selection valid detection events to a simultaneous detection of m 
photons at 5i (or Tti) and one photon at 82 (or 7 ^ 2 ), whereas simultaneous detection events of 
m photons at 5i and m photons at n\ or 1 photon at 52 and 1 photon at 7^2 are not taken into 
account. The corresponding setup is displayed in Eig. 2. Eor this configuration the {m+ l)-th 











Fig. 2 . Scheme of the modified setup. In contrast to Fig. 1 , m detectors i = 1 ,... ,m, 
are located in the far field at ri and one detector at r2. Assuming that each de¬ 

tector registers coincidentally exactly one photon the {m -f l)-th order correlation function 

^^TLS (■' 1 1 •'2) is measured. 


order correlation functions reads [19] 

(5i, ^2) (1 + ^2 cos(5i - Sz)) 

(TTi, TTz) (1 + 1^2 cos(5i - 52)) 

G^^tls (5i, ^ 2 ) (1 - 1^2 cos(5i - 52)) 

Gi^LS (^1 ,52) (1 - cos(5i - 52)) , 


where the visibility is given by ~ and 22 / is a constant given by 22 / = 

%t+i) (n)'”^*. Again, summing over all second order correlation functions yields 

the normalization factor 

L G'™+yO'G')=4^, (20) 

i=Si ,71 1 
j=^,7l2 

SO that the corresponding joint two-photon and single photon detection probabilities, derived 
from Eqs. (19) and (20), read 

Pf)„(5i,52) = i(l + r2(^+')cos(5i -52)) 

pf )„(7r,, 712) = ^(1 + r 2 ^rlP cos(5i - 52 )) 

, ( 21 ) 

pf )„(5i, 712) = - (1 - cos(5i - 52)) 

Pf )„(7i, , 52) = i(1 - cos(5i - 52)), 

and 

^2 rLs(5i) = P 2 TLs(^) ~ 2 ■ 

In this case the Bell inequalities Eq. (2) take the following form 
— 1< (^cos(5i — 52 ) — cos(5i — 52 )+cos(5[ — 52 )-Fcos(5( — 52 )^ — ^ < 0 . (23) 







It can easily be seen that, after replacing ^ Spe and assuming spe = this 

equation matches with Eq. (13) for two SPE. Moreover, a violation of Eq. (23) can be achieved 
if 'f' 2 'TLs'' ~ >1/This is realized for m > 5 in which case we have '^ 2 ^ 2^15 ~ 0-714 > 

1 /\/2. This violation of Bell’s inequalities demonstrates that measuring higher order correlation 
functions can indeed create and herald non-local multiparticle entanglement among photons. 
The interpretation for the considered setup is that a bunch of m > 5 photons at position 5i is 
entangled with one photon at ^ 2 , all m + 1 photons being emitted by classical light sources, in 
this case by two statistically independent TLS. 



Fig. 3. Experimental setup to measure the {m + l)-th order correlation function 
{81,82) with 2 Thermal Light Sources (TLS). For details see text. M: mirror; L: 

lens. 

The entanglement of the photons can be understood also from a different perspective. The 
detection of the first m photons can be interpreted as a projection of the two TLS sources onto 
a state which emits the subsequent photon not uniformly but in a strongly correlated manner, 
such that the Bell inequalities are violated. This projection of the sources can be seen as similar 
to Einstein Podolsky Rosen steering by photon subtraction as recently discussed in [20]. Prom 
this point of view the entanglement and violation of Bell’s inequalities is a consequence of 
the nonvanishing cross correlations between noncommuting quadrature phase components of 
the two spatially separated fields of the type (djdi)p(„,), with (dj)p(ra) = (di)p(m) = 0, where 
is the density matrix of the field after m photons have been recorded. Indeed, for the cross 
correlation coefficient = (alai){alao)p{m){d\di) [20] we obtain ^ = 

"^2TLs'^ ^ which approaches unity in the limit m 

4. Experimental setup and results 

To measure G^ 2 ^jll{d\^& 2 ) with two statistically independent TLS we used a mask with two 
identical slits of width a = 25 pm and separation d = 200pm. The mask is located a few 
centimeters behind a rotating ground glass disk which is illuminated by a linearly polarized 
frequency-doubled Nd;YAG laser at A = 532 nm (see Pig. 3). Within each slit a large num¬ 
ber of time-dependent speckles is produced, generated by the stochastically interfering waves 
scattered from the granular surface of the ground glass disk, representing many independent 
point-like sub-sources equivalent to an ordinary spatially incoherent thermal source. The co¬ 
herence time of the pseudothermal sources depends on the rotational speed of the disk [21] and 
was chosen in the experiment to Tc « 50 ms. To ensure that the two TLS radiate with equal 
intensity we enlarged the incident laser beam by a microscope objective to 1 cm providing a 























































Fig. 4. Visibilities normalized (m + l)-th order correlation func¬ 
tions (ri,r 2 ) = (ri,r 2 )/((Gyy^ 5 (ri))'”Gyy^ 5 (r 2 )) for two Thermal Light 

Sources (TLS) for m + 1 =2,...,9. The blue (dotted) curve shows the theoretical predic¬ 
tion ~ (solid) curve the experimental results. In the inset is shown the 

normalized correlation function ’ ’’ 2 ) m -f 1 = 7 as a function of the position of 

the detector at r 2 for ri = 0 (blue (dotted) curve: theory; red points: experimental results). 
The error bars reflect the statistical uncertainties. 


homogeneous illumination of the mask. To determine (5i, ^) we measured the intensity 

of the light scattered by the two TLS at different positions using a conventional digital camera 
placed in the focal point (Fourier plane) of a lens behind the mask (z « /) thus fulfilling the far 
field condition (see Fig. 3). Hereby, each pixel of the camera serves as a detector to register the 
intensity at position xjjz^ 9j, and we correlated m — 1 pixels at xi ~ 0i, each separated by one 
pixel along the y-direction to use m — 1 different pixels with identical xi-values, with another 
pixel at X 2 ~ 02- With more than a million of pixels the digital camera has the advantage that 
the amount of data accumulated in one frame to correlate the intensities at m different pixels 
is exceedingly larger than using m single photon detectors [22]. In order to obtain interference 
signals of high visibility, the integration time of the camera T, was chosen much shorter than 
the coherence time of the TLS, in our case T,- « 3 ms << T^. 

Fig. 4 displays the experimentally observed visibilities of the normalized {m+ l)-th order 
correlation functions g^'”y^y(ri,r 2 ) = G^"yyy(ri,r 2 )/((Gyy^ 5 (ri))'"Gyy^ 5 (r 2 )) for m-|-1 = 
2,... ,9. For illustration, the measured normalized correlation function ^ func¬ 

tion of the position of detector D„,+i at r 2 (for ri = 0) is shown in the inset of Fig. 4 and 
compared to the theoretical prediction of Eq. (19). It can be seen that the experimental results 
are in excellent agreement with the theory, in particular that the experimentally obtained vis¬ 
ibilities follow the expression ~ However, one can also see that experimentally 

for m = 5, due to remaining experimental uncertainties, a violation of Bell’s inequalities is not 
yet obtained, in contrast to the theoretical prediction = 0.714 > l/v/2. Nonethe¬ 
less, a visibility of = 0.743 ± 0.027 is achieved experimentally for m = 6. According to 

Eq. (23) this visibility demonstrates a violation of the Bell inequalities. 



















5. Conclusion 


In conclusion we showed that spatially separated photons emitted by statistically independent 
classical sources may display spatial correlations strong enough to violate Bell’s inequalities. In 
contrast to SPE [12], the violation does not occur for two recorded photons but for m>5 pho¬ 
tons registered in one mode and one photon in another mode. Due to experimental uncertainties 
the violation was demonstrated experimentally for m > 6. The spatial correlations among the 
photons can be understood from the perspective of state projection where the detection of the 
first m photons projects the sources onto a state which emits the subsequent photon in a strongly 
correlated manner such that Bell’s inequalities are violated. From this perspective the observed 
entanglement and violation of Bell’s inequalities is due to cross correlations among noncom¬ 
muting quadrature phase components of the two spatially separated fields only appearing after 
m photons have been recorded. 
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